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We introduce a charge coupled device (CCD) camera based detection scheme in 
dynamic light scattering that provides information on the single-scattered auto- 
correlation function even for fairly turbid samples. It is based on the single 
focused laser beam geometry combined with the selective cross correlation 
analysis of the scattered light intensity. Using a CCD camera as a multispeckle 
detector we show how spatial correlations in the intensity pattern can be linked 
to single and multiple scattering processes. Multiple scattering suppression 
is then achieved by an efficient cross correlation algorithm working in real 
time with a temporal resolution down to 0.2 seconds. Our approach allows 
access to the extensive range of systems that show low-order scattering by 
selective detection of the singly scattered light. Model experiments on slowly 
relaxing suspensions of titanium dioxide in glycerol were carried out to establish 
the validity range of our approach. Successful application of the method is 
demonstrated up to a scattering coefficient of more than fi s = 5 cm -1 for the 
sample size of L — 1 cm. 
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1. INTRODUCTION 

Dynamic Light Scattering (DLS) analyzes the intensity 
fluctuations of light scattered from a medium in the weak 
scattering limit. This is typically done by means of the 
normalized intensity autocorrelation function (ACF): 



(I(q,t)I(q,t + T)) t 



(i) 



where (•••)* denotes time averaging, r is the lag time and 
q is the scattering wave number or momentum transfer 
defined by scattering angle 9, the wavelength in vacuum 
A and the solvent refractive index n: q = Ai:nsin{6/2)/\. 
The measurable quantity g 2 (q, t) can be linked to actual 
physical microscopic properties by the normalized field 
autocorrelation function gi(q, r) via the Siegert relation: 



g 2 {q, t) = l + (3\gi{q,T)\ 2 , 



(2) 



where the coefficient j3 depends on the detection optics. 

Quite generally the the held autocorrelation function 
<?i(9> t) provides access to thermally driven local dynamic 
properties on length scales of the order 1 / q. A prominent 
example is the Brownian motion of colloidal particles in 
a solvent such as water. For this most simple case the 
normalized field correlation function can be written asi 
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gi(q,r) = exp(-D a q 2 T) = exp(—r/T c ), 



(3) 



where r c = 1/Dq 2 is the relaxation time and Dq is 
the particle diffusion coefficient defined by the Stokes- 
Einstein relation: 
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with i] the solvent viscosity, T is the sample temperature 
and R the particle radius. The equation JHJl is widely 
used in dynamic light scattering for the sizing of small 
particles. 

An essential condition for traditional dynamic light 
scattering to work is absence of multiple scattering. As 
soon as higher order scattering becomes considerable 
(typically if transmission in line-of-sight is lower than 
95 per cent) the measured intensity correlation function 
starts to deviate from the theoretical expectations, lead- 
ing to a faster decay. Furthermore, information on the 
scattering wave number as well as on the scattering angle 
6 is lost since the detected light is composed of several 
scattering events with unknown momentum transfer. 

The influence of multiple scattering can be reduced by 
decreasing the concentration of the sample under study 
or the cell size or by refractive-index matching 2 . The 
latter one is usually not possible without changing other 
sample properties. Limitations to the size of the con- 
tainer are set by the optical quality of the sample cells 
and by boundary effects. For cylindrical cells minimal di- 
ameters of typically 1 — 3 mm are used whereas in flat or 
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rectangular containers even smaller photon path lengths 
are achieved^. To maximally reduce the photons path 
lengths one can use fiber optical probes^ directly im- 
mersed in the (liquid) sample. This approach known as 
Fiber Optical Quasi Elastic Light Scattering (FOQELS) 
has been applied in a number of recent studies (see, e.g. 
RefsiWi?). The application of FOQELS is however lim- 
ited to backscattering angles around 180° and moreover 
the interpretation of the data is often complicated due to 
the incomplete suppression of multiple scattering. 

Quite a different way of actively dealing with multi- 
ple scattering has been put forward over the last two 
decades. The idea is to carry out two simultaneous 
DLS experiments with exactly the same scattering vec- 
tors in the same scattering volume and analyze the time 
cross-correlation function. It has been clearly shown 
that under proper conditions (see Refsi^iiSiil) the cross- 
correlation function equals the auto-correlation function 
for single scattering within the range of experimental 
resolution. Successful implementations of this scheme 
have been reported by several groups. The techniques 
are called two color DLS (TCDLS)ii*i£ and three dimen- 
sional DLS (3DDLS)2ii£ respectively and the latter one 
is available commercially^. 

Another cross-correlation approach uses a single-beam 
two-detector configuration Suppression of multiple 

scattering is based on a consequence of the van Cittert- 
Zernike theorem^, which states that intensity correla- 
tions in an observation region are closely related to the 
Fourier transform of the intensity distribution across the 
source. This means that a small region of single scat- 
tering (e.g. the volume of a focused beam) will produce 
large correlated areas (speckles), whereas a comparably 
large halo of multiple scatte red photons will give rise to 
small speckles (see Fig. 1(a) L This is reflected in the well 
known expression 16 for the speckle size 6x in the far-field 
geometry: 5x — Xz/S, where z is the distance from the 
light emitting object to the detector and S is the lateral 
extension of the object along one chosen direction. The 
consequences for a light scattering cross-correlation ex- 
periment are obvious. Spatially resolved detection of the 
scattered intensity will carry selective information about 
the spatial distribution of light in the scattering volume. 

If we assume that the size of single-scattering volume 
Si is equal to a average radius of beam cross-section 
of S\ ~ 20 um and the scattering mean free path be 
I w 1 mm then the volume of double-scattering extends 
over roughly a 25 times larger cross section. In prac- 
tice the dimension of the detected scattering volume will 
determine S2 both for weak and moderate multiple scat- 
tering. 

Scattered intensities measured in two points separated 
by distance Ax 
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will thus be corr elated only due to single scattering as 
shown in the Fig. 1(b) so the normalized intensity cross 
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Fig. 1. (a) Illustration of the van Cittert-Zernike the- 
orem: small coherently illuminated areas (such as a fo- 
cused laser beam) produces large correlated areas (speck- 
les) and and vice versa: large coherence areas (halo from 
multiple scattering) produces small speckles, (b) Sug- 
gested suppression scheme using cross-correlation pro- 
cessing: I) intensity values measured within the same 
speckle are correlated II) two different speckles are uncor- 
related III) for a superposition of small and large speckles 
intensities detected at a certain distance will be corre- 
lated only due to the larger speckles. 



correlation function (CCF) 

(I(q,t,0)I(q,t + T, Ax)) t 



(I(q,t,0)) t {I(q,t,Ax)) t 



(6) 



will provide the proper estimate of auto correlation func- 
tion of singly scattered intensity. Such an approach has 
already been successfully demonstrated by Meyer et aJj-L 
with the scheme based on cross-correlation of scattered 
intensities detected by two spatially separated fibers. 
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While the underlying optical background (van Cittert - 
Zernike theorem) is highly plausible it is more compli- 
cated to put forward a detailed theoretical description 
since this requires modeling of the low order scattering 
processes. Such treatment has been derived by Loclsii 
for the case of double scattering. He finds the multiple 
scattering suppression ratio to be approximately propor- 
tional to the speckle size ration 5*2 / Si if the detectors are 
placed at the distance Ax « \- z/S\. The suppression 
ratio can be improved by choosing a larger separation Ax 
albeit at the cost of a decreased signal. Choosing a large 
distance on the other hand might prove unnecessary for 
small amounts of multiple scattering. It is due to these 
practical difficulties, that the technically simpler single- 
beam cross-correlation geometry is often considered in- 
ferior to the two-beam realization (where a sample inde- 
pendent accurate theoretical description is available). 

Here we propose an extension of the single-beam cross- 
correlation method that allows to overcome this short- 
coming. Using a charge coupled device (CCD) camera as 
a detector we can analyze speckle correlations and adapt 
Ax thus assuring single-scattering detection with high 
accuracy. As we will show this flexibility together with 
the intrinsically high statistical accuracy of multi-speckle 
detection leads to a much improved performance of the 
single-beam two-detector configuration while essentially 
preserving its technical simplicity. 

2. EXPERIMENTAL SETUP 
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Fig. 2. Experimental setup. Light scattered at an angle 
of 90° inside the sample cell of thickness L = 1 cm passes 
a diaphragm and is collected by a CCD digital camera. 
In opposite direction light is collected by a mono-mode 
fiber and a photon counting module to be processed by a 
hardware correlator. The transmitted collimated inten- 
sity is measured to determine the scattering parameters 
Ms- 



The fluctuations of scattered light intensity are moni- 
tored in a traditional scattering geometry using a fixed 
scattering angle. As a light source solid-state laser system 
(Verdi Coherent Inc., USA) operating at the wavelength 
A = 532 nm is used. The beam is strongly focused by a 
lens with focal length / = 50 mm inside the sample cell 
which results in an average beam radius of 20 pm inside 
the scattering volume. The focal spot was positioned in 
the center of a rectangular quartz cell (Hellma GmbH, 
Germany) with inner base dimensions of 10 x 10 mm 
and height 45 mm. The light scattered from the sample is 
recorded simultaneously by single-mode fiber connected 
to a photon counting module (PerkinElmer, Canada) and 
analyzed with a digital correlator (Correlator.com, USA) 
from one side and a CCD camera from the opposite side. 
The CCD grayscale camera "Pixelfly" produced by PCO 
Computer Optics GmbH, Germany is configured to oper- 
ate in VGA mode (640 x 480 pixels) with 12 bits resolu- 
tion and 50 frames per second speed with exposure time 
To = 20 ms. This camera is based on a Sony ICX414AL 
image sensor with the pixel size of 9.9 x 9.9 pm. The 
dynamic range of analog/digital conversion of the CCD 
signal is 68.7 dB (according to manufacturer specifica- 
tions). A diaphragm with an opening diameter of roughly 
3.5 mm selects the range of scattering vectors seen by the 
CCD matrix and also determines the effective scattering 
volume. The distance for the beam center to the CCD 
matrix was z = 11 cm which corresponds to the speckle 



size of ca. 14 pm or 1.4 pixels. 

Scattering angles covered by the CCD chip are found 
in the range = 90° ±3.15°. The corresponding scatter- 
ing vector is equal to q = (24.59 ± 0.67) ■ 10 6 m -1 with a 
deviation of the order of 2.75 % from the average 90° an- 
gle. An improved angular accuracy could be achieved by 
placing the camera further away from the sample how- 
ever at the cost of decreased statistical accuracy due to 
the smaller number of detected speckles. A better way 
to deal with this problem would be the processing of an 
angle resolved digital image. However this has not been 
realized in this study. 

In parallel the intensity of the collimated beam trans- 
mitted through the cell is measured by a laser power- 
meter FieldMax (Coherent Inc., USA). The temperature 
is monitored by a digital thermometer with the probe 
placed close to the cell surface. 

As a model system we studied the Brownian motion of 
commercial Ti02 particles in pure (99.5%) glycerol. So- 
lutions were sequentially passed through filters with pore 
diameters 5 pm and 1.2 pm, respectively. Measurements 
on samples highly diluted with water reveal a mean hy- 
drodynamic diameter for TiC-2 particles of 293 ± 15 pm, 
in qualitative agreement with electron micrographs. For 
all measurements the temperature was in the range 20 ± 
0.5 °C. Corresponding variations of solvent viscosity are 
of the order of ± 5 We select Glycerol as a solvent 
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to decrease the particle diffusion coefficient Dq and thus 
enable real-time detection and processing of the scattered 
intensity fluctuations with a CCD camera 1 . 

Typically f 4 - fO 5 frames were collected correspond- 
ing to a measurement time of 1 5 — 30 minutes. The ac- 
tual number of collected correlation coefficients depends 
on the processing scheme (see the next section for de- 
tails) but usually is of the order of 10 s for the correlation 
coefficient of smallest delay time. 

The turbidity of the sample was characterized by 
means of the scattering coefficient p s . It is proportional 
to the particle density p and the scattering cross section 
a s : Us — pas ■ Due to the unknown amount of particles 
lost in the filtering process however no accurate density 
values are available. Since even fairly small particle den- 
sities lead to considerable multiple scattering we can nev- 
ertheless safely assume to work in the highly dilute limit. 
Experimentally p s can be estimated from the transmit- 
ted intensity It according to Lambert-Beer's law for the 
case of non-absorbing particles: It/h = Aexp(—p s L), 
where Iq is the intensity incident on the cell, L is the cell 
thickness, and A describes loss and deflection of intensity 
at the surface of the cell. The latter is independent of the 
particle density and was estimated from the transmission 
coefficient of a cell containing pure glycerol. 

3. PROCESSING TECHNIQUES 

A key element of our study is the optimization of multi- 
speckle detection and processing schemes. Our goal is 
to combine intelligent optical realizations with the power 
of modern parallel processing of a large amount of data 
acquired by a digital camera. Such optimized data anal- 
ysis can be both achieved in the time and and the space 
domain. As the first step we will discuss the applica- 
tion of the multi-tau scheme as developed by Schatzel to 
our multi-speckle analysis (previous realization for digi- 
tal cameras are described e.gi^). Secondly, we will show 
that a conceptually very similar approach can be used 
in the space domain in order to improve both the mul- 
tiple scattering suppression efficiency and the processing 
speed of our CCD detection scheme. 

A. Multi-tau correlation scheme 

The multi-tau correlation scheme was originally proposed 
by Schatzel^ to increase the accuracy of hardware corre- 
lators for large lag times. In most cases the linear spac- 
ing of lag times is not required to analyze single scatter- 
ing correlation functions. Instead one can increase the 
distance between lag points in the correlation function 



1 We note that Glycerol is hygroscopic and it is thus difficult to 
know precisely the exact water content. As a consequence a 
small uncertainty remains with respect to the solvent viscosity. 
We have assured, however, that all samples were prepared under 
identical conditions. Thus any systematic shift in the solvent 
viscosity will not affect the results of our study. 
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Fig. 3. (a) Principles of the multi-tau correlation scheme. 
Sequential intensity values I(t) are integrated for the 
computation of larger lag times, (b) Simplified object 
scheme to illustrate the cascading of linear correlators 
used for the realization of the multi-tau correlation 
scheme. 



for large lag times, thus saving valuable processing time. 
Such an approach can be easily and efficiently realized 
with a sequential doubling of the effective exposure time, 
which furthermore improves statistical accuracy. 

Let us suppose that the camera provides data with 
an initial exposure time equal to To with 1/tq Hz fre- 
quency which implies the absence of delays between the 
collection of sequential images. Then the doubling of the 
exposure time to 2to can be realized by integrating two 
sequential values of data obtained with exposure time tq. 
In the same way the effective exposure time can be in- 

From 



creased to 4tq and so on as it is shown in Fig. 3(a) 
the obtained time series of intensity fluctuations with dif- 
ferent exposure times the correlation coefficients can be 
calculated with a simple linear scale processing. 

In terms of software realization this is an ideal case 
for an object-oriented approach. The main object in 
this scheme is an Elementary Linear Correlator shown in 
Fig. |3(b)[ which receives the sequential data point 
as an input at each cycle i, multiplies the new value with 
previous data points I(i — k) for different linearly spaced 
integer delays k = 0, 1,2, ... and updates the mean val- 
ues of the corresponding products (I(i)I(i — k)). Ev- 
ery second cycle the correlator estimates the mean in- 
tegrated value of intensity for two sequential time steps 
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[I(2n) + I(2n+l)]/2 (where n = 0,1,2,...) and sends it to 
the next linear correlator, thus forming a cascading line 
of correlators. Since the data rate on a sequential corre- 
lator is half of the data rate compared to the preceding 
one, the evaluation period doubles at every correlator. 
Extra effort has to be taken to optimize the cycles with 
respect to limitations in computational power. 

The elementary correlator unit used in this study was 
designed to calculate the cross-correlation products to- 
gether with the auto-correlation. ACF coefficients are 
calculated on every pixel individually to obtain local val- 
ues of (I(x,t)I{x,t + r)) t and (I(x,t)) t . For the CCF 
with separation Ax we obtain (I(x, t)I{x + Ax, t + r)) t 
and (I(x,t)) t and (I(x + Ax, t)) t ■ These values are accu- 
mulated in each correlator for further spatial averaging. 

The scheme described here allows us to evaluate the 
data in real time and therefore does not set any limita- 
tions on the duration of the measurements. 

B. Binning technique 

The low dynamic range of CCD cameras in compari- 
son to photodiodes or photo-multipliers used in previous 
studiesii*^* is the main challenge for their utilization in 
a cross-correlation scheme which requires resolving the 
small signal from single scattering hidden by a dominat- 
ing signal from multiple scattering. The problem can 
be partially solved with an original binning technique 
developed in the course of this study. If the size of a 
speckle from single scattering exceeds the area of sev- 
eral pixels it can be approximated with integral values of 
these pixels intensities. Let us call the bin or meta-pixel 
the area S x x S y pixels represented by a single intensity 
value obtained by integration (or floating-point averag- 
ing) of intensities of included pixels. Due to multiple 
sampling the noise of measurements will be reduced by 
factor 'SxSijQ and thus the dynamic range defined as 
DR = 201og 10 SNR dB, where SNR is a signal-to- noise 
ratio, will increase on 10 log 10 S x S y dB. For example for 
a 10 x 2 window the dynamic range of our camera model 
is increased up to 81.7 dB or a SNR = 1.22 * 10 4 . The 
loss of spatial resolution reduces the statistical accuracy 
and intercept /3 only if the binning area S x x Sy is 
comparable or larger than the coherence area \J S\ x 52. 
As a consequence binning leads to partial suppression of 
multiple scattering by averaging out small speckles. Bin- 
ning introduces a two-dimensional filtering with a box- 
like kernel function which will reduce more efficiently the 
high-frequency spatial fluctuations (multiple scattering 
speckles) than the lower-frequency fluctuations that are 
connected to single-scattering speckles. In other words 
the binning technique in our cross-correlation approach 
can be considered a spatial analogue of the multi-tau 
technique in the time domain described above. 



very efficiently used to increase the statistical accuracy 
of a measurement by averaging the data along the bins. 
For the case of independent speckles different pixels or 
bins can be treated as separate photodetectors and thus 
their measurements can be processed all together to es- 
timate the correlation function. Actually two possible 
ways of dealing with this data exist. As it was men- 
tioned above correlators collect the time-averaged prod- 
ucts {I{x,t)I(x,t + r)) t as well as the mean intensity 
(I(x,t)) t for certain bins. Thus the further averaging 
followed by a normalization with the mean intensity can 
be performed with products and mean intensities as was 
proposed in the original multi-speckle scheme 22-23 : 



9 a 2 d (r) 



((I(x, t)I{x, t + r)) t ) x «/(£, t)I(x, t + r)) t ) 



where (. . .) x denotes the averaging along the entire two- 

2 

dimensional CCD matrix and I = ((I(x,t)I(x,t + 
T ))t)x — ((I(x,t)I(x,t + r)) x ) t is a mean value of inten- 
sity averaged over time and space. This "average-and- 
divide" sequence can be applied for the ideal case of a 
uniform illumination of the CCD matrix, i.e. when the 
mean intensity (I(x,t))t is not a function of x. Because 
for t — > oo, one finds I(x,i),I{x,t + r) — > (I(x,t))t and 
(I(x,t)I(x,t + r)) t — > {I(x,t))1, so for this processing 
scheme 



lim g 2 (r) 



f 



([(I(x,t)) t -If 
f 



-1 = K+l, 
(8) 

where K = ([I(x, t) — I] 2 ) x /I is the contrast of the pic- 
ture that would be obtained with infinite exposure time. 
For the case of uniform illumination K = due to the 
fact that (I(x, t)) t = I all along the detector matrix and 
thus §2 (t — ► oo) will approach 1 so the normalized field 
autocorrelation function <? 2 (t) — g 2 (t) — 1 will decay to 
zero. But even for a slightly non-uniform illumination 
when K > the non-zero additive component appears in 
the measured values g\ (r) . 

To overcome this problem the original multi-speckle 
scheme was somewhat modified in our study by normal- 
ization of the locally estimated correlation coefficients, 
and averaging only after that ( "divide-and-average" se- 
quence) : 



92 a (r) 



(I(x,t)I(x,t+T)) t 

COM))? 



(9) 



Still this requires the deviation of the mean intensity 
along the matrix to be small in comparison to the noise 
level. For a perfectly uniform illumination of the matrix 
both algorithms provide the same result. 



C. Multi-speckle averaging technique 

The ability of cameras to register a large number of in- 
dependently fluctuating speckles simultaneously can be 



4. RESULTS AND DISCUSSION 

We have carried out a series of experiments using samples 
of different scattering strength, hence different amounts 
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Fig. 4. Left: Representative area taken from the 
recorded images for a weakly scattering (A) and a 
strongly scattering (B) sample. Right: Spatial intensity 
fluctuations along the vertical columns of the CCD ma- 
trix. Sample A: /j, s = 0.101 cm' 1 , estimated speckle size 
Sx = 28.02 pixels, sample B: fi s = 5.92 cm -1 , Sx = 0.76 
pixels 



of multiple scattering, ranging from the dilute limit 
fi s L = 0.1 to a regime /i s L — 5.92 where the transmitted 
beam is attenuated to only 0.27% of its initial intensity. 

A. Spatial intensity correlations in the speckle pattern 

As soon as multiple scattering effects become consider- 
able the detected intensity distribution will consist of a 
distribution of correlation lengths. The longest corre- 
lation length, Sx — \z w 300 um ~ 30 pixels can be 
associated with single scattering whereas higher order 
scattering leads to smaller speckles. Figure 0] illustrates 
this by showing the intensity fluctuations along the CCD 
matrix columns. For the case of single scattering sam- 
ple (/x s = 0.1 cm , upper plot) the anisotropy of the 
speckle pattern can be clearly seen. The reason for the 
anisotropy is the horizontal confinement of the incident 
focused beam that defines the scattering volume in this 
regime. The scale of intensity fluctuations along the "x" 
dimension are of the order of tens of pixels so the cor- 
related areas of intensity or speckles are large. With in- 
creasing scattering coefficient fi s the fluctuations become 
more pronounced and the correlation length decreases. 
The horizontal speckle-size is always given by the dimen- 
sions of the detected scattering volume both in the single 
and multiple scattering case. 

For a quantitative estimation of the actual speckle size 
the normalized spatial auto-correlation function can be 
used: 



ff2(Aa;) 



(I{x)I(x + Ax)) Xjt 
(I(x)) Xtt (I(x + Ax)) x y 



(10) 



where (• • • ) Xj t denotes the time and space averaging. In 
practice the averaging is performed through all pixel pairs 
separated vertically by a distance Ax. A total of 8 
independent speckle images is analyzed with an expo- 
sure time of 20 ms. Temporal fluctuations are much 
slower, typically of the order of 0.3 seconds or more, 
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Fig. 5. Spatial auto-correlation for different sample tur- 
bidities as a function of vertical pixel separation Ax: (o) 
H s = 0.10 cm- 1 , (■) fi„ = 0.39 cm" 1 , (y) Ma = 0.88 
cm" 1 , (♦) (M s = 1.43 cm- 1 , (A) ju s = 2.54 cm" 1 , (*) ^ s = 
5.92 cm -1 . Dashed line ( ) represents the same esti- 
mated in horizontal dimension. Inset shows speckle sizes 
Sx (o) and Sy (■) as a function of /i s for both vertical 
and horizontal separations respectively. The maximum 
amplitude of the correlation function is limited to about 
0.5 since the speckle size in y-direction is comparable to 
the pixel size 



and could thus be neglected. Estimated values for the 
mean speckle size Sx, Sy can be defined by the condition 
g s 2 (Ax, Ay<$:Sy)-l = 1/e, g s 2 (Ax < Sx, Ay) - 1 = 1/e. 
Since speckles in vertical direction are always small, the 
limit Ay -C Sy cannot be reached in our setup and we 
have thus calculated the correlation length in Ax direc- 
tion using the correlation function normalized with re- 
spect to g%(Ax = I, Ay = 1). Spatial ACFs for samples 
of various turbidities are shown in Fig. [S] together with 
the estimated speckle sizes. In vertical direction, with the 
increase of fi s , the speckle size decreases rapidly. In the 
multiple scattering limit the speckles are almost isotropic 
with a size of approximately lxl pixels which means 
that the initially focused beam is dispersed into a diffuse 
scattering cloud. 

The spatial resolution in our experiments allows direct 
estimation the relative weight of single scattering contri- 
butions. For separations much larger than 3 pixels the 
cross correlation signal will be dominated by single scat- 
tering. A separation comparable to the size of a single 
scattering speckle, Ax ~ 30 pixels is a good compromise 
value separating single from multiple scattering. 

Note that using a CCD camera we are however not 
restricted to fixed detector positions as in previous stud- 
ies. Smaller separations Ax could be used in the weak 
scattering limit whereas a gradual increase of Ax in the 
strong scattering limit simultaneously optimizes multiple 
scattering suppression and the signal-to-noise ratio. 
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Fig. 6. Auto correlation functions obtained with dif- 
ferent binning areas together with the cross-correlation 
function. As the vertical size of the binning area in- 
creases the ACF approaches the CCF indicating partial 
but not sufficient suppression of multiple scattering . The 
inset shows the estimated particle size (o) as a function 
of binning size. The actual size is indicated with by a 
dashed line. 



B. Pixel binning 

Instead of pixel by pixel cross correlation we can apply 
the superior binning approach described before. Proper 
selection of the binning areas can increase the signal-to- 
noise ratio and also partially filter out the small multiple 
scattering speckle. The latter effect is demonstrated in 
Fig. El for different numbers of binning pixels in the verti- 
cal direction. Since the horizontal speckle size is constant 
the binning in this dimension was selected to be 2 pixels. 

Auto correlation functions calculated from such larger 
areas are found increasingly closer to the single scatter- 
ing function as the binning area is increased. However 
complete suppression is difficult to achieve since multiple 
scattering suppression scales linearly with the number of 
pixels whereas cross correlation scales exponentially with 
pixel separation. Another disadvantage of very large bin- 
ning areas is the reduction of the statistical ensemble. 
Ideally the binning area is chosen somewhat smaller than 
the size of single scattering speckle in order to retain a 
high number of independent speckles. In our case the 
best choice in x-direction is in the range of typically 5-20 
pixels as it can be seen from Fig. [5] 

C. Intensity Correlation Functions 

Fig. shows the correlation function for one of the most 
strongly scattering samples using different measurement 
schemes. With our optimized cross-correlation (Ax = 40 
pixels) and binning approach (20 x 2) the measured corre- 
lation function perfectly agrees with the single scattering 
result. 

A set of measurements has been carried out and the 



Fig. 7. Correlation functions obtained by different means 
for an essentially multiple-scattering sample (/i s =3.93 
cm -1 ). The ACF obtained from the hardware correla- 
tor ( ), the ACF from CCD detection with no bin- 
ning (□) and the CCF with a separation Ax = 40 pixels 
(o) are shown together with the ACF obtained for singly 
scattering sample (•). Inset shows the intercepts for ACF 
(□) and CCF with separation Ax — 40 (•) pixels as func- 
tions of /i s . 
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Fig. 8. Particle size D determined by different methods 
as a function of sample turbidity. The CCF diameter re- 
mains unchanged up to the limit strong multiple scatter- 
ing. T. ACF from the hardware correlator (•) and from 
ACF without binning (□) show a rapid decrease as the 
scattering coefficient /x s increases. The shaded area indi- 
cates the experimental value expected for single scatter- 
ing (including the uncertainty in particle size and solvent 
viscosity) 
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measured particles diameters are presented in Fig. [5] as 
a function of scattering coefficient fj, s . For the case of 
dilute samples {^i < 0.5 cm -1 ) all processing schemes 
yield the same results within the experimental error. In- 
creasing fi s results in smaller apparent values of the par- 
ticle size for the case of autocorrelation measurements 
( ACF) . At the same time the particle diameter obtained 
using our multiple scattering suppression scheme (CCF 
processing) stays the same for all n s < 5.5 cm -1 corre- 
sponding to fiL — 5.5. Only for the most turbid sample 
fj, s ps 5.92 cm -1 a noticeable difference is observed due to 
the loss of the single scattering signal and the unavoid- 
able increase of error. 

Finally we would like to comment on the performance 
of the multi-speckle averaging scheme implemented in our 
approach. Neglecting any angular dependence we are 
treating all x-columns equally. As we have seen the size 
of a single scattering speckle roughly amounts to 30 x 2 
pixels. Since the full CCD chip has 640 x 480 pixels we 
detect approximately 5000 independent speckles. This 
number can be increased only if the diameter of the in- 
cident beam is increased, but at the expense of a inferior 
performance in multiple scattering suppression. 

5. CONCLUSION AND OUTLOOK 

CCD-camera based light scattering can be used to effi- 
ciently suppress the undesirable effects of multiple scat- 
tering. Our single-beam cross correlation scheme com- 
bines for the first time the advantages of multiple- 
scattering suppression and multi-speckle detection. The 
latter reduces the collection time for slowly evolving sam- 
ples dramatically and provides the necessary statistical 
accuracy to detect even small cross correlation signals 
deep in the multiple scattering regime. Successful ap- 
plication of these combined processing schemes has been 
demonstrated for sample turbidities as large as fi s L ps 5. 
We think that our relatively simple approach, heavily re- 
lying on a flexible processing scheme, can be very useful 
for a variety of applications for the study of complex flu- 
ids and soft materials. With the maximum frame rate 
available in our experiments the most appropriate areas 
of application can be found in the field of slowly relax- 
ing systems such as glasses and gels or dense surfactant 
solutions, to name a few. But at the current rate of devel- 
opment of image sensors and computer hardware much 
smaller lag times seem to be feasible in a near future. 
Already now, with the use of modern fast complemen- 
tary metal-oxide semiconductor cameras with embedded 
programmable digital signal processors the real-time on- 
camera calculation of the cross-correlation functions is 
theoretically possible which in turn would tremendously 
improve the time resolution. 
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List of Figure Captions 

Fig. 1. (a) Illustration of the van Cittert-Zernike the- 
orem: small coherently illuminated areas (such as a fo- 
cused laser beam) produces large correlated areas (speck- 
les) and and vice versa: large coherence areas (halo from 
multiple scattering) produces small speckles, (b) Sug- 
gested suppression scheme using cross-correlation pro- 
cessing: I) intensity values measured within the same 
speckle are correlated II) two different speckles are uncor- 
rected III) for a superposition of small and large speckles 
intensities detected at a certain distance will be corre- 
lated only due to the larger speckles. 
Fig. 2. Experimental setup. Light scattered at an angle 
of 90° inside the sample cell of thickness L = 1 cm passes 
a diaphragm and is collected by a CCD digital camera. 
In opposite direction light is collected by a mono-mode 
fiber and a photon counting module to be processed by a 
hardware correlator. The transmitted collimated inten- 
sity is measured to determine the scattering parameters 
Ms- 

Fig. 3. (a) Principles of the multi-tau correlation 
scheme. Sequential intensity values I(t) are integrated 
for the computation of larger lag times, (b) Simplified ob- 
ject scheme to illustrate the cascading of linear correla- 
tors used for the realization of the multi-tau correlation 
scheme. 

Fig. 4. Left: Representative area taken from the 
recorded images for a weakly scattering (A) and a 
strongly scattering (B) sample. Right: Spatial intensity 
fluctuations along the vertical columns of the CCD ma- 
trix. Sample A: /j, s — 0.101 cm -1 , estimated speckle size 
5x = 28.02 pixels, sample B: fi s — 5.92 cm -1 , Sx = 0.76 
pixels 

Fig. 5. Spatial auto-correlation for different sample tur- 
bidities as a function of vertical pixel separation Ax: (o) 
H s = 0.10 cm" 1 , (■) n s = 0.39 cm" 1 , (y) Us = 0.88 
cm" 1 , (♦) (i s = 1.43 cm" 1 , (A) [i s = 2.54 cm" 1 , (*) ^ s = 
5.92 cm" 1 . Dashed line ( ) represents the same esti- 
mated in horizontal dimension. Inset shows speckle sizes 
5x (o) and 5y (■) as a function of /j, s for both vertical 
and horizontal separations respectively. The maximum 
amplitude of the correlation function is limited to about 
0.5 since the speckle size in y-direction is comparable to 
the pixel size 

Fig. 6. Auto correlation functions obtained with dif- 
ferent binning areas together with the cross-correlation 
function. As the vertical size of the binning area in- 
creases the ACF approaches the CCF indicating partial 
but not sufficient suppression of multiple scattering . The 
inset shows the estimated particle size (o) as a function 
of binning size. The actual size is indicated with by a 
dashed line. 

Fig. 7. Correlation functions obtained by different means 
for an essentially multiple-scattering sample (fi s =3.93 
cm" 1 ). The ACF obtained from the hardware correla- 
tor ( ), the ACF from CCD detection with no bin- 
ning (□) and the CCF with a separation Ax = 40 pixels 
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(o) are shown together with the ACF obtained for singly 
scattering sample (•). Inset shows the intercepts for ACF 
(□) and CCF with separation Ax = 40 (•) 
Fig. 8 Particle size D determined by different methods 
as a function of sample turbidity. The CCF diameter re- 
mains unchanged up to the limit strong multiple scatter- 
ing. T. ACF from the hardware correlator (•) and from 



ACF without binning (□) show a rapid decrease as the 
scattering coefficient /i s increases. The shaded area indi- 
cates the experimental value expected for single scatter- 
ing (including the uncertainty in particle size and solvent 
viscosity) 



